Real-space renormalization yields finite correlations 



Thomas Barthel,' Martin Kliesch,' and Jens Eisert'-^ 

'institute for Physics and Astronomy, Potsdam University, 14476 Potsdam, Germany 
^Institute for Advanced Study Berlin, 14193 Berlin, Germany 
(Dated: March 01, 2010) 

Real-space renormalization approaches for quantum lattice systems generate certain hierarchical classes of 
states that are subsumed by the multi-scale entanglement renormalization ansatz (MERA). It is shown that, 
with the exception of one spatial dimension, MERA states are actually finitely correlated states, i.e., projected 
entangled pair states (PEPS) with a bond dimension independent of the system size. Hence, real-space renor- 
malization generates states which can be encoded with local effective degrees of freedom, and MERA states 
form an efficiently contractible class of PEPS that obey the area law for the entanglement entropy. It is shown 
further that there exist other efficiently contractible schemes violating the area law. 
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I. INTRODUCTION 

Renormalization group (RG) methods aim at solving many- 
body problems by treating energy scales in an iterative fash- 
ion, progressing from high to low energies |[l[ 01- One of 
its earliest formulations is the real-space RG which works by 
repeated steps of thinning out local degrees of freedom and 
rescalingof the system as in Kadanoff's block spin transfor- 
mation yh. In real-space RG approaches to quantum lattice 
models fllHl, in each RG step r, the system is partitioned 
into small blocks. From those blocks high-energy states are 
eliminated and the Hamiltonian Hr+i for the renormalized 
system is obtained by applying the coiTesponding projection 
operators, exactly Hr+i = -FV+i^^t-PJ+i or in some appro- 
priate approximation, followed by a coarse-graining of the lat- 
tice. This is iterated, e.g., until a step t = T is reached where 
the renormalized system consists of a single small block for 
which the ground-state |gSj.) can be obtained exactly. Apply- 
ing the RG transformations in reverse order yields an approx- 
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P^lgSrp) to the ground-state of the original 
model. Those states, generated by the real-space RG, fall into 
the class of so-called tree tensor networks (TTN) A re- 
cent more elaborate real-space RG scheme, the multi-scale en- 
tanglement renormalization ansatz (MERA) ||7[[^, a genuine 
simulation technique for strongly correlated systems, allows 
in each RG step for local unitary operations to be applied be- 
fore the elimination of block basis states lf36ll . The technique 
generates hence a more general class of (variational) states, 
referred to as MERA states; see Fig.[T] 

Whereas the degrees of freedom of MERA and TTN states 
are organized in a hierarchical structure encoding correlations 
on different length scales, there exists a different class of so 
called ^«;refy correlated states where the degrees of freedom 
are organized in a strictly local manner. For D = 1 dimen- 
sional systems they are often referred to as matrix product 
states 1191 4 1 ill , and for Z? > 1 as tensor product ansdtze or 
projected entangled pair states (PEPS) |[l2l - [l6ll : see Fig. |2] 
PEPS are the basis of powerful numerical techniques, such as 
the extraordinarily successful density -matrix renormalization 
group method iflTlfTsll . 

In this work, we establish the surprising fact that, for 
D > 1, real-space RG, despite of the inherently hierarchical 
nature of the procedure, generates states that capture corre- 



lations by local degrees of freedom. More specifically, it is 
shown that MERA states form a subclass of PEPS, unifying 
both approaches. This also explains the failure of real-space 
RG for some situations for which merely anecdotal evidence 
had previously accumulated. 

PEPS, TTN, and MERA are all tensor network states 
(TNS). In terms of an orthonormal product basis |cr) = 
^i=i Wi) foi" ^ lattice of N sites, TNS are of the form 
J2(T V'ct|c) where the expansion coefficients tAo- are encoded 
as a partially contracted set of tensors; Fig. |2] Recently, this 
notion has been generalized to the fermionic case lfl9l422ll . 
For a PEPS, to each site i, a tensor Ai is assigned which has 
one physical index Ui and further auxiliary indices - one for 
each nearest neighbor - which need to be contracted to obtain 
ijjcr'. Fig. 12] For TTN and MERA, the tensors are arranged in 
a hierarchical pattern with the physical indices in the lowest 
layer; Fig.[T] The number of degrees of freedom of a TNS can 
be tuned by changing the number x of values each auxiliary 
index runs over. Increasing x for a fixed structure of the TNS, 
enlarges the variational space, allowing for a more precise ap- 
proximation to the exact ground-state in a variational method, 
but increases computation costs. Hence, x is called the re- 
finement parameter of the TNS. The computational costs for 
efficient simulation techniques scale polynomially in x- 

In this article the following is shown. 

• For D > 1 spatial dimensions, MERA states with a 
refinement parameter x can be mapped efficiently to 
PEPS such that the resulting PEPS refinement param- 
eter xpeps is a system-size independent function of x- 

• For _D = 1, TTN and MERA states can in general not 
be encoded efficiently as PEPS. For a fixed x there are 
MERA states with an entanglement entropy that scales 
logarithmically in the system size, as occurring in criti- 
cal models ll23ll24ll . In this sense MERA states are more 
useful than PEPS for this case. 

• For D > 1, there exist other efficiently contractible 
TNS, based on quantum cellular automata, that gen- 
erate an amount of entanglement violating the area law. 

The ability to map MERA states efficiently to PEPS for £> > 1 
implies that D > 1 MERA states always obey the entangle- 
ment area /aw just as PEPS ll23l - t25ll . This behavior is shared 
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Figure 1: A ID MERA with linear branching ratio b — 2. Circles, 
squares, and the triangle denote tensors, the lines denote contractions 
of those tensors. The squares are (partial) isometrics that map two lo- 
cal subsystems HI and "HJ+i into one ''^ Kadanoff 's block 
spin transformation. The circles denote unitary operators, disentan- 
glers, that reduce the entanglement between T-LJ ® "HJ+i and the 
rest of the system before the action of the isometry. When, for the 
mapping to a ID PEPS, MERA tensors are assigned to lattice sites 
according to Eq. ([6}, as in the diagram, stacks of tensors occur: There 
are sites where the number of assigned MERA tensors diverges with 
the system size, implying that the mapping is inefficient. The stack- 
ing of tensors can be avoided by choosing tensor positions according 
to Eq. l|7j. 

by ground-states of non-critical systems and critical bosons. 
Ground-states of critical fermions, however, can violate the 
area law ll23i l24l l26l - {30ll . Consequently x needs to be scaled 
polynomially in the system size in order to describe such crit- 
ical fermionic systems accurately. Otherwise, the real-space 
RG schemes addressed here 10, IH I?!] are imprecise in that 
case. The remaining advantage of D > 1 MERA is that local 
observables and correlation functions can be evaluated effi- 
ciently, whereas, for PEPS, approximations are necessary. In 
this sense, MERA states form an efficiently contractible sub- 
class of PEPS. 



II. GENERAL PROCEDURE FOR MAPPING TNS TO 
PEPS 

All TNS can be mapped to PEPS, although not necessarily 
in an efficient manner. To map a TNS to a PEPS one can 

1. assign each tensor of the TNS to a specific site of the 
physical lattice iIjtIi 

Vphys:={0,...,L-l}^cZ^, (1) 

and, 

2. for each contraction line that connects the tensors, de- 
cide on a specific path for that line on the edges fphys of 
the physical lattice, 

fphys {(r, r') G Vphys x Vphysl |r - r'|i = 1}, (2) 

see Fig.|2] The tensors composing the PEPS are then obtained 
by introducing for each edge of the lattice an auxiliary vector 
space that is the tensor product of the vector spaces of all TNS 
contraction lines that traverse that edge. The elements of the 




Figure 2: (a) Procedure for mapping a TNS (left) to a 2D PEPS 
(right), by assigning tensors to lattice sites and contraction lines to 
paths on the lattice, (b) The elements of the PEPS tensors are deter- 
mined by the elements of the tensors composing the TNS. 



PEPS tensor for site i are determined by the elements of all 
the TNS tensors that were assigned to site i. See Fig.|2j3. 

Applying this procedure for a ID MERA state inevitably 
results in a PEPS refinement parameter xpeps that diverges 
with the system size. This is not just a feature of the specific 
procedure. In Sec.[Vl a family of ID TTN states is constructed 
for which any mapping to PEPS necessarily requires xpeps to 
diverge with the system size. 



in. EFFICIENTLY MAPPING MERA TO PEPS FOR D>1 

Given a family of MERA states for different lattice sizes 
L, a mapping of the MERA states to PEPS is called efficient 
if there exists an upper bound xpeps on the resulting PEPS 
refinement parameter that is independent of L. 



A. Qualitative argument for the existence of an efficient 
mapping of D > 1 MERA to PEPS 

The following argument motivates why an efficient map- 
ping of MERA to PEPS should be possible for Z) > 1 but 
not for D = 1. Let us assign to each contraction line of the 
MERA state a finite cross-section, e.g., equal to a^~^ with 
the lattice spacing a. Then one can ask what ZJ-dimensional 
volume V{t) the contraction lines of a certain layer r con- 
necting to layers with t' < t cover. Those contraction lines 
of layer r have length £(r) oc alf , where b is the linear 
branching ratio of the MERA. The number of lattice cells in 
layer r is b^^~'^^^; Fig.[T] Hence, the volume covered by the 
contraction lines of layer t is V{t) cx a^^^l{T)b'^'^~'^'^^ oc 
^Z5T-(d-i)t rj-jjg density of the MERA contraction lines, or 
more precisely, a resulting estimate for the average number 
of contraction line paths traversing a unit cell of the physical 
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lattice (t = 0) is hence 



T=Q T=0 

't foi-L> = l 



(3) 



Note that for an edge that is traversed by n paths, one obtains 
an upper bound Xpeps = x" to the corresponding PEPS re- 
finement parameter, i.e., n = log^(xpEPs)- As T = Xog^^L, 
ID MERA with a fixed refinement parameter x have accord- 
ing to Eq. (O the potential to encode states with a logarithmic 
scaling of the entanglement entropy, as occurring in critical 
ID systems. See Sec.[V]for an example. For D > 1, however, 
Eq. (O suggests that there is enough space on the physical 
lattice to assign the MERA contraction lines to paths on the 
lattice in such a way that, for a fixed x, the resulting PEPS 
has a bond dimension Xpeps that is independent of the system 
size. That this is indeed possible is proven constructively in 
the following. 



B. Preconditions for MERA states 

In order to show that the mapping presented in the follow- 
ing is efficient, it is necessary to exploit the defining proper- 
ties of MERA states that correspond directly to features of the 
real-space RG and can be summarized as follows. 

1 . The MERA state is a TNS for a D-dimensional square 
lattice (Vphys, fphys) consisting of unit cells with 
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(4) 



2. The MERA consists of T layers of tensors labeled by 

T = i,...,r. 

3. There is an upper bound x on the dimension of the vec- 
tor spaces associated to the tensor indices, and an upper 
bound Co on the order of each tensor. 

4. With each layer, we associate a coarse-grained square 
lattice Cr of (L/fe^)^ cells of the physical lattice 



Cr ■■= {0,...,L/fe^- 1}^ C Z^, 



(5) 



and £o '■= Vphys. Every cell of lattice Cr contains cor- 
responding cells of lattice Cr-i Issll . 

5. There exists an assignment of the tensors of layer r to 
cells of the lattice Cr such that 

(a) the number of tensors inside a single cell is 
bounded from above by a constant Ct, 

(b) the distance of contracted tensors is bounded from 
above by Cr, where the distance of a tensor of 
layer r to a tensor of layer r' < r is defined as 
the Li distance of their corresponding cells in Cr 



6. For |r — r'l > Ct, there are no contractions between 
tensors of layer r with tensors of layer r' . 

The upper bounds x, Co, Ct, Cr, and Ct are required to be 
independent of the system size L. ll40ll The stated conditions 
guarantee that the MERA states feature a so-called causal 
cone i^. Hence, local observables can be evaluated efficiently 
if all tensors are chosen isometric. For the mapping described 
in the following, this is however not necessary. As we require 
only upper bounds on the MERA refinement parameter, which 
may hence actually vary from tensor to tensor, the apparent 
restriction to square lattices is not essential. The conditions 
stated above are met for all typical MERA structures consid- 
ered in the literature so far In Fig.|4^ a 2D MERA with & = 2, 
Co = 8, Ct ^ 2, and Ct = 1 is depicted. Fig.|4]3 shows a 2D 
MERA with = 3, Co = 10, Ct = 4, and Ct = 1. One can 
reach Cr = 2 for both examples. 



C. The efficient mapping 

Let us explain a general scheme for mapping MERA states 
for _D-dimensional systems efficiently to PEPS. The precon- 
ditions of Sec, nil Bl are assumed to be given. A simple proce- 
dure to assign the MERA tensors to certain lattice sites is to 
put the tensors of cell n E Cr of layer t to the site 



rr{n) = ITn e Vphys 



(6) 



as in Fig.[T| The problem with this approach is that one gen- 
erates stacks of tensors at certain lattice sites, i.e., there ex- 
ist positions r S Vphys to which a number of tensors is as- 
signed that is not independent of the lattice size. For ex- 
ample, at site r = (0, . . . , 0) a number of cx T = log^ L 
tensors accumulate. Further stacks of tensors with height 
oc T' accumulate for example at lattice sites with coordinates 
6"'" (1, . . . , 1) <E Vphys, see Fig.[T] It is necessary to avoid such 
stacks of tensors, because they imply in general that xpeps 
diverges with the system size. Stacks can be avoided by shift- 
ing the allowed tensor positions for different layers relative to 
each other. One possible such choice for rr{n) G Vphys is 



rr{n) = b'^n + b'^ with n £ Cr 



(7) 



and e := (1,...,1) G Z^. With this choice, tensors of layers 
r and t' can end up at the same site only if 

rr{n) = rr'{n') 

^ = b^'-^n'^ + fe^'-'-i - b-^ Vi, 

i.e., only if the two tensors belong to the same layer t' ^ t 
and the same lattice cell n' = n within that layer, following 
from the fact that the right hand side of the second equation 
is not an integer otherwise. The possible tensor positions of 
layers r form disjoint sublattices Vr of the physical lattice 

Vphys ■ 



Vr := {rr{n)\n G Cr} C Vphys, 

V^nVr'=0 ^r^r'- 



(8) 
(9) 
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Figure 3: Disjoint sublattices Vt C Vphys to which MERA tensors are 
assigned and disjoint subsets of edges £t C f^phys to which MERA 
contraction lines are assigned for b — 2. In our construction, the 
paths assigned to contraction lines from tensors of layer 2 to tensors 
in layer 3 are for example restricted to edges from £2 U £3. 



All coordinates r; of r € Vr have a b-adic valuation of r — 1, 
where the 6-adic valuation Vb{n) of an integer n is defined 
such that Vb{n) ~ r iff r is the largest integer such that n mod 
6^ = 0, for example, ■^2(12) = 2. 

Avoiding stacks of tensors is not sufficient for an efficient 
PEPS encoding. In D = 1, all contraction lines are assigned 
to paths that necessarily stack up on the x-axis. Fig. [T] This 
stacking of the paths can be avoided in Z3 > 1 by assigning 
contraction lines between tensors of layers r and r' to paths 
that are restricted to edges from certain subgrids £r and £t-' 
and that are shortest paths with respect to the Li distance on 
Vr U Vr'. Here, a grid £r is defined as the subset of phys- 
ical edges connecting nearest neighbors of the lattice Vr on 
straight lines; see Fig. [3] 

£r {(r,r + e,) G £phys| Vb{rj) = r - 1 V,^J (10) 

^ £rr\£r' =fl}yr^r' (H) 



[ei]j = Si J. 



with 

For this choice for tensor positions and paths of MERA 
contraction lines an upper bound for the resulting PEPS re- 
finement parameter xpeps follows: Contraction lines assigned 
to an edge e = (''", t') G £t contract tensors of layer t with 
tensors of layers r' where |t' — t| < Ct- For a layer r' with 
r' > T, tensors from at most {2Cr)^ cells of Cr' around the 
cell corresponding to site r can have contraction line paths 
traversing edge e. From the layers t' with r' < r, tensors 
of at most (2a)^ EfJo ^ells can contribute. Thus, the 
number of contraction line paths traversing edge e and hence 
log^ (xpeps) are bounded from above by 



log^(XPEPs) < (2a)'' (Ct + 6^(^-+i))aa. (12) 

As this upper bound is independent of the system size, the 
presented mapping of MERA to PEPS is efficient. 
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Figure 4: (a) Unit cell of a specific 2D MERA state with linear 
branching ration b — 2. With each layer, corresponding to a single 
renormalization step, unitary disentanglers are applied that reduce 
the entanglement between blocks of 2 x 2 sites with the rest of the 
system. Then, an isometry maps from those 2x2 sites (dots) into 
one (crosses), (b) Unit cell of an alternative 2D MERA state. In each 
layer, unitary disentanglers are applied that reduce the entanglement 
between blocks of 3 x 3 sites with the rest of the system. Then, an 
isometry maps from those 3x3 sites into one. (c) Mapping of the 
MERA state from (a) to a refined PEPS with 5t = 1. The diagram 
shows the assignment of two layers of the MERA, composed of dis- 
entanglers u and isometrics w, to the physical lattice. Each edge of 
the lattice is traversed by at most two contraction line paths. The 
resulting PEPS has xpeps = • 



D. Lower bond dimension by PEPS refinement 

As it stands, for each layer, all tensors of a given lattice cell 
n ^ Cr of layer t are assigned to the same physical lattice site 
r^(n) G Vphys according to Eq. (|2|. Therefore, a considerable 
number of contraction lines that start at the tensors of a given 
cell n may traverse the same edges around r^(n) and cause 
hence a relatively high xpeps- While this is unproblematic for 
the purpose of proving the existence of an efficient mapping, 
the situation can be improved for numerical purposes, e.g., by 
introducing for each site of the physical layer b^'^^ — 1 aux- 
iliary sites with 5t > 0, resulting in refined lattices Vpfjy,, and 
V'r- The corresponding refined PEPS has tensors for the phys- 
ical sites and tensors for the auxiliary sites, where the latter 
ones carry no physical indices. The sites from Vpjj^^ allowed 
for tensors of layer t are then defined as 



Trin, m) := b'^+^^n + b^m + h'^-^e 



(13) 



with n G £t and G {0, . . . , 6'^'" - 1}'' 



Lattice cells in layer t are again labeled by n G £r and m 
labels now the possible positions for tensors inside that cell. 

Fig. |4}; displays a mapping of the MERA from Fig. |4^ to 
a refined PEPS with 6t = \. Each edge of the grid is tra- 
versed by at most two contraction line paths. The result- 
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Figure 5: The ID TTN state discussed in Sec.[V]for T = log L — 5. 
For the chosen subsystems At, the entanglement entropy is logarith- 
mic in the system size, Sat = {T + 1) /2. 



ing PEPS has consequently Xp^ps = X ■ Due to the re- 
finement of the physical lattice, the PEPS consists however 
of l"l^pt,ysl &*^-°|Vphys| instead of |Vphys| tensors. A re- 
fined PEPS is transformed to a "normal" PEPS by contract- 
ing the PEPS tensors for the b^'^^ — 1 auxiliary sites with 
the tensor for the corresponding physical site, resulting in 

XPEPS — IXpEPS ) 



IV. EXAMPLES FOR PARTICULAR 2D MERA STATES 

The scheme displayed in Fig.|4}; for mapping the 2D MERA 
defined in Fig.|4^ to a PEPS results in the PEPS refinement pa- 
rameter XpEPs'^ = if one uses a refined PEPS with (5t = 1 
where isometrics are located at m = (1,1) and disentanglers 
at m = (0, 0), according to Eq. (fT3] l. Using, instead, the ten- 
sor coordinates according to the most simple scheme, Eq. (|7J, 
yields xpeps = X^- 

An analogous scheme for mapping the slightly more com- 
plicated MERA defined in Fig. to a PEPS results in the 
PEPS refinement parameter Xpeps'' ~ if one uses a re- 
fined PEPS with St = 1 where the isometrics are located at 
m = (1, 1), the 2 X 2 disentanglers at m = (0, 0), the 2 x 1 
disentanglers at m = (1,0), and the 1x2 disentanglers at 
m = (0, 1), according to Eq. (fTSl l. 



V. ID TTN AND MERA STATES CANNOT BE MAPPED 
EFFICIENTLY TO ID PEPS 

To show that a ID MERA can in general not be mapped to a 
ID PEPS with a bond dimension xpeps that is independent of 
the system size L, we construct a family of graph TTN states 
ifl lsill for which the entanglement entropy grows logarithmi- 
cally with L for a suitable bipartition of the system; see, e.g., 
Refs. LSjil^lllt] for a numerical analysis. We choose a TTN 
state, so a MERA state without disentanglers, with isometrics 
mapping from two qubits to one, i.e., x = 2 and b = 2. Each 
representative of the family of states is defined on L = 2^ 
sites {0, . . . ,L — 1} with a positive odd integer T. For sim- 
plicity, the TTN is embedded into the entire lattice of L sites. 



and each isometry is considered as a unitary having one input 
from the previous layer and one input |0); see Fig. |5] The 
state of the the top layer t = T and the two-site gates u are 
defined as 



|0) 



and u 



(14) 



X and Z denoting Pauli matrices, and |0) being an eigenstate 
of Z. For layers r = 1, . . . , T the state IV'r-i) is generated 
from |-0T-) by applying gates u to sites 2'^(fc — 1/2) — 1 and 
2'^fc — 1 for fc = 1, . . . , 2^^^. The entanglement is computed 
for the subsystem At consisting of the last px sites, where 
pi := 1 and pt+2 4pT — 1. Since all the gates, spec- 
ified in Eq. ( fT4l i. are mutually commuting, all gates that are 
supported entirely on At or entirely on its complement Aj^ 
can be disregarded for the computation of the entanglement 
entropy 



Tip At ^OgPAr with pAr = Tr_^J, |'0t)('0' 



The locations pr of the bipartitions are chosen such that, for 
odd T, exactly (T + l)/2 gates act across the cut. Since each 
such gate generates one pair of maximally entangled qubits, 
one obtains Sat = (T^ + l)/2 which is logarithmically diver- 
gent in the system size L and implies that any ID PEPS en- 
coding of the given graph state requires a Xpeps that diverges 
with L. 



VI. EFFICIENTLY CONTRACTIBLE TNS THAT VIOLATE 
THE AREA LAW 

As shown here, unlike for D = 1 spatial dimensions, 
MERA states for D > 1 obey the entanglement area law 
and not a log-area law, Sal = ^(L^~^ log-^), as it oc- 
curs for critical fermionic models 11231 l24l 1261 - BOll . This raises 
the question of whether any efficiently contractible tensor net- 
work automatically yields an area law. This is however not the 
case. In order to show this, we construct, for a _D-dimensional 
cubic lattice of sites, a family of efficiently contractible 
TNS based on a unitary quantum cellular automata (QCA). 
For a specific choice of the constituting tensors, one obtains 
instances that violate the area law for generic bipartitions of 
the system. 

Let us consider a QCA consisting of T layers t = 1, . . . , T, 
where each layer consists of two sublayers. With the first 
sublayer, Ki is applied which consists of (L/2)^ local uni- 
tary gates s supported on plaquettes of 2 x • • • x 2 sites 
each. The operator K2 for the second sublayer is identical 
to A'l except for a relative shifting of the gate positions by 
(1, . . . , 1). Therefore, periodic boundary conditions are im- 
posed, and L is required to be even. The initial state j^/'o) is 
a product state of plaquette states |0) for 2 x • • • x 2 

sites each, where the plaquette positions coincide with those 
of the gates in K2- The plaquette states are product states 
of 2^~^ maximally entangled pairs of qubits sitting each at 
the ends of the plaquette diagonals, e.g., \(j>) ~ uo.ilO)*^^ for 
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Figure 6: The graph state QCA discussed in Sec. IVII for D = 1 
and D — 2. Each layer of the QCA moves the maximally entan- 
gled qubits that initially are located on next nearest neighbor sites 
two steps further apart from each other. This increases the entangle- 
ment entropy for a given bipartition of the system in every step. For 
an appropriate choice of the number of layers, the states violate the 
entanglement area law while still being efficiently contractible. 



D = 1, where u is chosen according to Eq. ( fT4l i and the in- 
dices label the sites the gate acts on. For D ~ 2, the pla- 
quette state is 10) = U(o,o),(i,i)"(i,o),(o,i) lO)*^"*- The plaque- 
tte operators s, composing the Ki, are chosen as products of 
swap operators Sij\(Ti(jj) = IcjTj) that act similarly on the 
qubits at the ends of the plaquette diagonals. For example, 

s = 'S'(o,o),(i,i)5'(i,o),(o,i) for 13 = 2. 

The QCA layers and the initial state are invariant under 
translations by two sites and rotations by 7r/4 and so are all 

states \tpr) ■— {K2Kiy\ipo)- Like \i/jo), every state \^r) is 
a product state of 2^^^{L/2)^ maximally entangled qubit 
pairs. If two entangled qubits have positions r ± Ar in \ipr) 
there is exactly one corresponding entangled qubit pair at po- 
sitions r ± (1 + j^^)Ar in |?/v+i)- Applying one QCA layer 
after another, distances of entangled qubits increase by Ay/D 

in each step, e.g., -ftrii^2W(2,2).(3,3) = W(o,o),(5,5), see Fig.|6] 

For generic choices for bipartitions of the system into two 
parts, Al C Vphys and its complement, where yl^ is connected 
and has a volume cx L^, the corresponding entanglement en- 
tropy will violate the area law if the number of layers, T, is 
chosen appropriately. Consider as an example the bipartition 
with Al = {0,...,L/2 - 1} X {0,...,L - The 
subsystem boundary is formed by the planes {r\ri = 0} and 
{r\ri = L/2}. Each plane is crossed by a number of different 
pairs of entangled qubits that is proportional to its area and to 
T, as long as i/2 > 2T. Consequently, 



(15) 



For a choice T oc log L, this yields a log-area law Sal = 
il{L^~^ log L). But an upper bound on the computation 
cost for the evaluation of local observables with respect to 
QCA states of the given class (with arbitrary u) is of order 

0(22^" T), i.e., Cl(L2(i°si)""' logL) for T cx logL. This 



cost is not polynomial in L and the QCA are for this T hence 
not efficiently contractible in an obvious fashion. However, 
the computation cost is, of order 0{L^ {log Ly^^) for the 
choice T = (logi)^/^, i.e., polynomial in L. The resulting 
entanglement entropy is = fl{L^^^ (log Ly^^) which 
violates the area law by the sub logarithmic factor (log L)^/^. 
Note also that even a QCA consisting of a single layer of 
k X ■ ■ ■ X k plaquettes, where k is allowed to grow as k cx 
(log L)^/^, can also violate the entanglement area law, albeit 
only for specific bipartitions of the system, while being effi- 
ciently contractible. 



VII. CONCLUSION 

In this work, we have shown that MERA states for D > 1 
can be encoded efficiently as PEPS. From the perspective of 
numerical simulations for strongly correlated many-body sys- 
tems, this means that D > 1 MERA states form a subclass 
of efficiently contractible PEPS. From a physical perspective, 
the result implies that real-space renormalization techniques, 
despite the scale-invariant features of the TNS they generate, 
give rise to states that can be encoded with local degrees of 
freedom. As a corollary, it follows that D > 1 MERA states 
feature an area law for the entanglement entropy ll23|j25ll . 
Consequently, the refinement parameter x needs to be scaled 
polynomially in the system size in order to describe critical 
fermionic systems accurately. Otherwise, the real-space RG 
schemes addressed here H El 01 are imprecise for such sys- 
tems. Constructing explicit examples, it has been shown fur- 
ther that ID MERA states can in general not be encoded ef- 
ficiently as PEPS and that, for D > 1, there exist efficiently 
contractible TNS, based on quantum cellular automata, that 
violate the entanglement area law. It is the hope that this work 
further clarifies the intricate structure of the small number of 
effective degrees of freedom that typical ground-states of local 
quantum many-body systems explore. 
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